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ABSTRACT 

It is argued that in hot gauge field theories, "Hard Thermal Loops" leading order calcu- 
lations call for a definite sequence of angular averages and discontinuity (or Imaginary part 
prescription) operations, and run otherwise into incorrect results. The ten years old coUinear 
singularity problem of hot QCD, provides a striking illustration of that fate. 
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1. Introduction 

The intrinsic non perturbative nature of non zero temperature Quantum field Theories 
has long been recognized [1]. Naive thermal perturbation theory can nevertheless be devised, 
both in imaginary and real time formalisms [2], but then, it promptly appears that, under 
certain circumstances, the original perturbative series must be re-organized. Such an example 
of re-organization is provided by the so called Resummation Program [3]. This program, RP 
for short, is a resummation scheme of the leading order thermal fluctuations which, in the 
literature, are known under the spell of Hard Thermal Loops. Whenever one is calculating a 
physical process related to thermal Green's functions whose external/internal legs are soft, it is 
mandatory to trade the naive thermal perturbation theory for the RP. The softness alluded to 
above, refers to momenta on the order of the soft scale gT, where T, the temperature, stands for 
the hard scale and g for any relevant (bare/renormalized) and small enough coupling constant. 

The RP which has been set up in order to remedy an obvious lack of completeness of 
the naive thermal perturbation theory, has produced interesting results, but has also met 
serious obstructions in the infrared regime of the theories [4], [5]. Within the Resummation 
Program itself, the solutions proposed so far [5,6,7], however interesting in their own respect, 
cannot be organized a systematic way and display a crucial dependence on the process under 
consideration. 

In the past few years, another resummation scheme of the leading order thermal fluctua- 
tions has been proposed by the present author, and seems to avoid all of the infrared problems 
met by the standard RP, [8] . In loop calculations (where both resummation schemes are being 
used as effective perturbation theories, ruling the leading thermal fluctuations) the former only 
differs the standard RP, the sequence along which are performed an integration over a looping 
energy po, on the one hand, and a sum over the number N of HTL-seU energy insertions, on 
the other hand. This new resummation scheme, has accordingly been denoted PR for short, 
and is, by construction, a Perturbative Resummation scheme of the leading order thermal 
fluctuations. 

Such is not the case of the RP. While HTL n-point vertices are purely perturbative 
objects, effective propagators are not, giving rise to pole residues and dispersion laws that 
cannot be obtained out of pure thermal perturbation theory. In the infrared limit, in particular, 
those effective propagators are known to give rise to series in the coupling constant that are 
no longer Taylor's but Laurent's series [8]. In any higher number of loops calculations, this 
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simple, still crucial difference, allows one to foresee easily, why, in the RP case and not in the 
PR one, some enhancement mechanisms may occur (and effectively will!) so as to make higher 
number of loops as important, and even bigger, than lower number of loops diagrams [9]. 

In this note, we take advantage of a comparison of the two RP and PR resummation 
schemes to point out some overlooked aspects concerning "the historical derivations" of the 
famous collinear singularity problem of hot QCD, [5]. As many of our further calculations have 
displayed ever since, these aspects are generic and extend beyond the historical derivations to 
be recalled shortly. 

We will be using the convention of upper case letters for quadrimomenta and lower case 
ones for their components, writing, for example P = (po^p). Our conventions for labelling 
internal and external momenta can be read off Fig.l. 

2. The collinear singularity problem of hot QCD 

This almost ten years old issue is the following. The soft real photon emission rate out 
of a Quark-Gluon Plasma in thermal equilibrium involves, in particular, the calculation of the 
quantity 



The discontinuity is to be taken in the energy variable po, by forming the difference of R and 
yl-indiced P-dependent quantities, and within standard notations, fermionic HTL self energies, 
effective propagators and vertices are respectively given by 




(2.1) 



*5,(P) = 



a = R,A 



(2.2) 



(P) + ieaPo 




(2.3) 



%(Pa,Q/3,P5) = -^e (7m + 




{Pa,Qp,P^)) 



(2.4) 



(2.5) 



where K is the hghtlike four vector (1, A;). In view of (2.4), four terms come about, three of 
them proportional to a coUinear singularity. These singular terms are the two terms with one 
bare vertex 7^, the other V^'^^ , plus the term including two HTL vertices, V^'^^ . Thanks 
to an abelian Ward identity peculiar to the high temperature limit, a partial cancellation of 
these collinear singularities occurs, but out of the term including two T^^^ vertices, a coUinear 
singularity remains. 



\J 47r Q.K + ie J J (27r)3 

X [Tr [*SAiP)^) - Tr (^*Sr{P')^)] (2.6) 

where, the soft photon being real, Q is the lightlike 4- vector Q= qQ =q{l, q). In the literature, 
this result is written in the form 



^ s=±l,V=P,P' 



C^* f d^P 
e 



where the overall 1/e comes from a dimensionally regularized evaluation of the factored out 
angular integration appearing in (2.6), and where iSgiV) is related to the effective fermionic 
propagator usual parametrization [2], 

*5„(P) = l:Y. fsAsipo + iea,p) (2.8) 

^ s=±l 

with Ps = {l,sp), the label s referring to the two dressed fermion propagating modes. One 
has then 

As{po + iea,p) = A"(po,p) = as{po,p) - ine{ea)l3s{po,p) (2.9) 
where e(x) is the distribution "sign of x". 



3. On an improper derivation 

The third term, with the two HTL vertices, T^^^, involves two angular integrals Wi that 
are, with i G {1, 2}, 

W.(P, P'^=fff^-f ^ , (3.1) 

' J Atv J 4n {K-P + i€)(K-P' + ie){K'-P + i€)(K'-P' + i€) 



While W2 only is met in the standard RP calculation, Wi and W2 come into play within the 
PR scheme. As shown in [10], these integrals are given by very cumbersome expressions, so 
that, following [5], we will illustrate our point on the two diagrams including one bare vertex 
7^, the other V^^^ . The R/A real time formalism conveys us to the expression 

f dk Tr (*Sr{PW*Sr{P')^) 

discp / ^ (3.2) 

J 47r {K-P + i€){K-P' + ie) 

where the superscript (★, 1) in the left hand side refers to a self energy diagram involving 
two effective propagators and one vertex HTL correction. The steps leading to the coUinear 
singularity of (3.2) are as follows (the second reference of [5], using the R/A real time formalism, 
is followed here). The discontinuity in po is taken in consistency with the R/A formalism, by 
writing 



Tr (*Sr{P)^*Sr{P')$:) Tt S r{P)$'' S r{P')$) 

discp — ^ = — ^ 

{K-P + ie){K-P' + ie) {K-P + ie){K-P' + ie) 

Tr (*Sa{.PWSr{P'W 



(3.3) 



{K-P -ie){K-P' + ie) 

Closing the po integration contour in the upper half complex po plane, one selects a pole term 
coming from the vertex T^^^ which reads 

/ if FFT^ {'SAPWS.iP'W) (3.4) 

Since P' = P + Q, we have indeed, as a building bolck of (3.4), the expression 

hAiPWSRiP'w) (3.5) 
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The angular average therefore develops a collinear singularity in a neighbourhood of K=Q, 
not compensated for by the numerator, and whose "residue" reads 



C(^*) 6{Q-P) 



Tr (^*SAiPWSRiP')^) (3.6) 



The first factor, singular at e = 0, is obtained by using a dimensional regularization of the 
angular integral. Up to regular contributions that we do not consider here, this singular piece 
of (3.5) translates, for H^'^'^^Q), into the singular result 

r(*.*;l).r.^ _ ^/(^^^eW f 

(27r) 



nr'^(Q) = -V-— / — ^(l-2r^^(po)) 



X diQ-P) Tr (^^SAiPWSniP')^) (3.7) 



The contribution of (3.7) to the soft photon emission rate being proportional to its imaginary 
part, the emission rate is thus plagued with a coUinear singularity. In its original (published) 
version [5], this result is obtained within the imaginary time formalism, continued to real 
energies, where the RP has been first devised [3]. 

In a recent article [10], the analogous PR calculation has been proven to be mass/coUinear 
singularity free. In particular, such is the case of any of the contributions to the soft real photon 

emission rate due to H^'''^ '^\Q) quantities, the sum of which, on and A^' corresponding 
to II^^'*'^\Q), as depicted on Fig.l. Thus, at face value, a collinear singularity shows up in a 
resummation scheme, not in the other. Let us now examine how can this be so. 

In a PR scheme, one has indeed for the potentially dangerous most part of any H^'"^ '^^ 
quantity, the expression [10], 



ieV{ .. ) J](m^)(^'+^-^"'--) / -|^(l-2nHPo)) 

n,n' 



^ \iP^ + iepo)'' J J 47r K-P + ie K-P' + ie ^ ' 

where the discrete sums on n and n' extend over a finite set of integers and need not be furhter 
specified here. Performing the angular integration, one finds 



J 4^ 



ln^!±|^^±M (3.9) 



47r K-P + ie K-P' + ie 2Q-P + ieq P^ + iepo 

The collinear domain is included in the phase space region where Q-P ~ 0, where one has 

1 . P^ + 2Q-P + iep'o ^ 1 1 2Q-P + ieq 

2Q-P + ieq^ P^ + iepo P^+iepo 2{P^+iepo)^ ^' ^ 

Obviously, this HTL-vertex induced behaviour is potentially mass singular by the light cone 
region P^ ~0. This means that in a neighbourhood of Q-P ~ 0, (3.9) mixes up with partial 
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propagators, S^^^ ^\P{P')), own potentially mass singular behaviours. Recalling that we 
have 

the whole integrand structure of (3.8), in the coUinear domain, is readily seen to "boil down" 
to a simple shift of inverse power 

^ ^ M = N + N' (3.12) 



(P2 + iepQ)M i^p2 + igpo)M+l ' 

As demonstrated in [10], similar shifts, as well as many more transformations, are proven to 
leave totally unaffected the robust mass singularity cancellation patterns which guarantee the 
regular character of the PR calculation. 

That is, in contradistinction with the RP scheme, a picture emerges out of the PR scheme, 
where effective vertices potentially mass singular behaviours, melt with partial effective propa- 
gators own potentially mass singular behaviours into structural patterns which rule the overall 
compensation of octwo/ mass/collinear singularities, [8,10]. 

However, it is important to remark that, apart from their constitutive difference, not 
exactly the same steps have been followed in either RP and PR schemes. Both resummation 
schemes are here being developed within the R/A real time formalism, so that, starting from 
the PR expression (3.8), exactly the same procedure as followed from (3.3) to (3.7) may be 
applied, with, as a straightforward result 



Q{st) „2_2 



em, 



( .. ) J](m^)(^'+^-^-'-^-) / -^(l-2nHPo)) ^ 
tt' J (27r)3 p'p 

'2v^2 t ulWri' r d2v2 



[ -p'^T^up'r [-p^si(p)] - 

(P'2 + iep'o)^' (P2 - iepoY 



In words, (3.8) is now discovered to display the same coUinear singularity as the one plaguing 
the RP result (3.7), in sharp contradiction with our claim that any H^'^ '^\Q) contributes 
mass/collinear singularity free quantities to the soft photon emission rate. 

The contradiction is of course only apparent and entirely rooted in the fact that passing 
from (3.1) to (3.3), the prescription of discontinuity inpo has been commuted with the angular 
integration on K. As we have just come to see, in both RP and PR schemes, the effect of that 
improper commutation (the discontinuity of the integral is prescribed by the formalisms, and 
not the integral of the discontinuity!) is to provide the coUinear singularity with a somewhat 



absolute status. In either schemes in efTect, the colhnear singularity just factors out, in total 
independence of the remaining integrations to be performed. 

The lack of commutativity of the discontinuity and angular integration can be the matter 
of a direct observation. From (3.9), one has 



47r K-P + ie K-P' + ie Q-P + ieq q P^ + iepo 

i7re(po)e(-P2) 



(3.14) 



Q-P + ie 

In the right hand side, the term proportional to the distribution S{Q-P), has coefficient 
zero, and not the colhnear singularity C^^^^/e, as historically derived whithin the reverse, 
improper sequence. Though the illustration above, (3.5)-(3.6), is given within the R/A real 
time formalism only, it must be noticed that, mutatis mutandis, both the historical derivations 
of [5] hinge on a reverse sequence calculation, 

/ ^ f{discp , or Im} ^^^) = -i {i ,or ]} Uk (3.15) 

7 47r V ^K-P + ie) K-P' + ie 2' J K-Q + ie 

along which a colhnear singularity is identified at K=Q. Then, as advertised after (2.7), one 
relies on a dimensionally regularized angular integration to find out the singular counterpart 

-i {\ .or-}! dk = -i^S{Q-P) (3.16) 

where e is the regularizing parameter of an angular integration performed at D = 3 + 2e spatial 
dimensions. This derivation, however, is inaccurate. This can be verified by a direct calculation 
of (3.15)'s right hand side, which gives instead 

f ^ fdiscp , ' = "9(-f;)^g-?> (3.17) 

7 47r V K-P + ieJ K-P' + ie Q-P-ie 

A first important remark is that under subsequent angular (p) and energy (po) integrations, 
the above reverse sequence result (3.17) will lead to more singular expressions than the proper 
sequence one, (3.14). Likewise, the correct prescription of +ie is not preserved by the sequence 
(3.17). Now, the point is that nothing like the "absolute", overall factoring out colhnear 
singularity of (3.16), i.e. the term C^^^-'/e, ever appears in the final result, and this underlines 
the illicit character of the manipulations leading to it. 

Discarding thus (3.16), one may still compare the right hand sides of (3.14) and (3.17). 
Then, apart from the fact that a non vanishing result is obtained for strictly colhnear p 
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and q momenta only, another essential difference comes about immediately, which is the sign 
distribution e(po)- In real time formalisms, sign distributions have long been noticed to be 
crucial so as to preserve integrable and non-integrable mass singularity compensations [11], 
and as observed here again, they are a natural outcome of the calculational operations proper 
sequence (3.14). 

Indeed, this is the very place where intricacies are taking place, between angular averages 
peculiar to n-points iJTL-vertices, and discontinuity or Imaginary part prescriptions. Staring 
at (3.9), for example, one can observe how the angular average is able to correctly reproduce 
all of the relevant internal legs RjA specifications, the itq, iepo and iepg, which do not appear 
at all in the left hand side. As can be read off [10], sections 4 and 5, this is not an isolated 
situation and rather reveals to be general a mechanism. 

That is, the proper sequence of angular average and discontinuity operations complies 
with the required, exact sign distributions which otherwise, are clearly endangered. 

By the way, this elucidates also the intriguing point made in the third reference of [8], 
after equation (3.32), which was then left as an issue. There, the discontinuity in po of the 
angular integral W2 was taken to be 



J 47r K-P' + teJ 



5{K-P) f dK' {K-K 



/\2 



47r {K'-P + i€){K'-P' + i€) 



whereas it is manifest that in conformity with the common use, the reverse improper sequence 
had been followed. In order to preserve mass singularity cancellations though, the correct 
sign distribution, eipo), had to be restored by hand (this restoration was motivated a heuristic 
way, by recalling that the vertex corrections F^-^^, were after all nothing but leading order 
approximations of full order 3^ expressions, endowed with explicit iepo specifications). 

Since further, higher number of loops calculations [12] enforce this peculiarity too, it seems 
reasonable to conclude that in any RP or PR calculations, and any real or imaginary time 
formalism, angular averages should definitely be performed first. 



4. The coUinear problem and the proper sequence 



A full, extensive treatment of the soft real photon emission rate RP calculation, as the 
proper sequence is followed, falls far beyond the scope of the short present note, and could be 
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postponed to a future publication. We will here content ourselves with an illustration of how 
the coUinear singularity problem gets translated when the angular integration is performed 
first. 

From (3.2) and (2.8), one obtains 
iImU^*^*''\Q) = -ieV [ ^{l-2npipo)) ^ E l3sipo,p)f3s'iPo,p) 



/■ d*P 1 

J i^T^r pp 

X E (s'(3Ap'o,p') (1 - s^) discp AUP) ln(^^) + (P ^ P')) 

s,s'=±l ^ P PO P J 

X (1 - .^) discp A|,(P) ^ In ^ (4.1) 

where the distribution /3s{po,p) of (2.9), though textbook material, may be worth recalling 
here in view of the forthcoming comments 

f3s{Po,p) = Zsip)S {po - UJs{p)) + Z_s{p)5 {po + UJ-s{p) ) 



m 



2 (i_5a.)e(-p2^ 



+ 1^ ^ ^-^ 2 (4-2) 

(p(l -sf)-^ ((1 - sf) In + 2.)) + ^(1 - sfr 

It is easy to check that the first line of (4.1), with the two functions (3s, is safe. For 
example, one may rely on the integration contour technique used for energy sum rules, [2], 
chapter 6, to show that po and p integrations are well defined. Likewise, inspection shows that 
no problem is inherited from the second and third lines of (4.1). At po = ±p, in effect, the 
potentially dangerous term of ln\po +p/po — p\ is, according to the value of s = ±1, either 
depleted by its own squared expression, or cancelled by a factor of (1 — spo/p)^- 

For the last term of (4.1), a potentially singular most contribution to Im II^^'*'^\Q) is 
generated by taking the discontinuity in po of the first factor, AJ^(P). One gets 



2 2 

-e m 



/-^(l-2n^(po)) E i^-s'^)PAPo,P') 

^"^^^ s,s'=±l P 



X i^-^j)f3siP0,p)^ In^ (4.3) 
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The Dirac pieces of both spectral densities /?s and /Jg/, pose no problem, whereas the other 
parts, proportional to the Heaviside distributions 0(— P^) and ©(— P'^), and hereafter denoted 
by Ps,s', yield 



2 2 

-e m 



5:(l-A^.(po,p)^ln^ (4.4) 



s=±l 



In the phase space domain where 2Q P is vanishing, the integrand of (4.4) behaves like 

(l-2n^(pz)) (i-s'^±^)^,,{q+pz,p')x ^ ^.(p^,p) (4.5) 

s'=±l P P s=±l ^^^^ 

where we have introduced the scaling variable z = po/p, and where p''^{z) = + 2qpz + q^. 
Writing 

1 - 2nF{pz) = 2[nF{p) - npipz)] + (1 - 2nF{p)) (4.6) 

the last term only, in the right hand side, can possibly induce a singularity by the light cone 
region (z = ±1), reached from below, and this would be singular contribution can be expressed 
as 



27T 



G(-P2) J pdp{l- 2nF{p)) dz ^ 



E (l + + VPZ,p'ivz)) (4.7) 



77=±1 

where use has been made of the relation (3s{—po,p) = P-s{po,p)- Now, against all odds, the 
light cone {1 — z ~ 0) behaviours of (3-{pz,p), (l + ^"'^T^ j , and /?_ {q + r]pz,p'{r]z)), make it 
straightforward to check that the above potentially dangerous most contribution is a regular 
one indeed. 

Since, out of the last term of (4.1), no other sign of a possible singularity shows up, it 
appears that the two n^'*'^-* (Q) contributions to the soft real photon emission rate are regular, 
in contradistinction with the usual, singular results. 
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5. Conclusion 

In this note, we have stressed that whatever the formaUsm (real or imaginary time) and 
the resummation scheme (standard (RP) or perturbative (PR)), a definite sequence of angu- 
lar averages and discontinuity (or Imaginary part) operations, must definitely be preserved. 
Moreover, this proper sequence is nothing exotic, as it simply corresponds to the sequence 
naturally prescribed by any real or imaginary time formalism one may use. 

Ignoring that proper sequence may lead (and has led) to incorrect derivations and/or 
results. In particular, infrared singularity cancellation patterns have long revealed to be par- 
ticularly sensitive to the angular averages and discontinuity operations proper sequence. This 
is also a point which, we think, should be kept in mind within the context of numerical simu- 
lations. 

For the hot QCD coUinear problem met in the RP scheme, this point certainly matters 
a lot. Our rapid analysis of the one effective vertex topologies, n^*'*'-*^), has not allowed us to 
detect any relic of a mass/collinear singularity in the soft photon emission rate, the proper 
sequence being followed. 

Now, as advertised in section 3, this indicates only that the proper sequence is less singular, 
under subsequent integrations, than the reverse improper one, and let totally open the issue 
concerning the two effective vertex topology, n*^*'*'^), contribution to the emission rate. 

In other words, the whole correct settings of the hot QCD coUinear singularity RP- 
problem.. if any! .. has to be worked out again. 
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Figure caption 



Fig.l: A graph denoted by {N,N'; 1), with N{N') insertions of HTL self energy along 
the P(P')-line, one bare vertex —iej^, and one HTL vertex correction (2.5). It is a PR scheme 
object. In the standard RP scheme, the two internal P(P')-lines are to be replaced with the 
full effective propagators *S'a(P(P')) of Equation (2.2), and this corresponds to the two RP 
graphs denoted by (★, 1) in the text. 
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